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A major obstacle toward a quantitative verification, by comparison to experiment, of the theory of 
the excitations of the fractional quantum Hall effect has been the lack of a proper understanding 
of disorder. We circumvent this problem by studying the neutral magneto-roton excitations, whose 
energy is expected to be largely insensitive to disorder. The calculated energies of the roton at 
1/3, 2/5 and 3/7 fillings of the lowest Landau level are in good agreement with those measured 
experimentally. 

PACS numbers:71.10.Pm,73.40.Hm 
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Quantitative tests of the theory of the fractional quan- 
tum Hall effect (FQHE) have focused in the past primar- 
ily on the gap to charged excitation, determined experi- 
mentally from the temperature dependence of the longi- 
tudinal resistance. A factor of two discrepancy between 
theory and experiment has persisted over the years, be- 
lieved to be caused by disorder for which a quantitatively 
reliable theoretical treatment is not available at the mo- 
ment. In recent years there has been tremendous experi- 
mental progress in the measurement of the energy of the 
neutral magneto-roton excitation M, both by inelastic 
Raman scattering Mm and by ballistic phonon absorp- 
tion |g-0|, and its energy has been determined at Lan- 
dau level fillings of 1/3, 2/5, and 3/7. While the neutral 
magneto-roton is of great interest in its own right, be- 
ing the lowest energy excitation of the FQHE state, the 
chief motivation of this work is the observation that the 
disorder is not likely to affect its energy significantly, in 
contrast to the energy of the charged excitation, because 
the roton has a much weaker dipolar coupling to disor- 
der due to its overall charge neutrality, and the coupling 
is further diminished because the disorder in modulation 
doped samples is typically smooth on the scale of the size 
(on the order a magnetic length) of the spatially localized 
roton. There is also compelling experimental evidence 
for the insensitivity of the roton energy to disorder: the 
same roton energy was found for samples for which the 
gaps in transport experiments differed by as much as a 
factor of two [p|. The roton therefore provides a won- 
derful opportunity for testing the quantitative validity 
of our understanding of the excitations of the fractional 
quantum Hall state. With this goal in mind, we have 
undertaken a comprehensive and realistic calculation of 
the roton energy at several filling factors in the lowest 
Landau level (LL). 

The neutral excitation of the FQHE will be treated 
in the framework of the composite fermion (CF) theory 
P-|lO|, the composite fermion being the bound state of 
an electron and an even number of flux quanta (a flux 



quantum is defined as (j)o — hc/e). According to this 
theory, the interacting electrons at Landau level filling 
factor v = n/{2pn ± 1), n and p being integers, trans- 
form into weakly interacting composite fermions at an 
effective filling ly* = n; the ground state corresponds to 
n filled CF-LLs and the neutral excitation to a particle- 
hole pair of composite fermions, called the CF exciton. 
Microscopic wave functions for the CF ground state and 
the CF exciton are readily constructed by analogy to the 
known wave functions of the electron ground state at fill- 
ing factor n, ^-^^ , and its exciton, <i>^^: 



$s 



$e 



--VlllY[{zj 

j<k 

Vlll W {zj - 
j<k 



Zk 






Zk 



)2P$; 



ex 

n 1 



(1) 



(2) 



where Zj — 



iyj is the position of the jih particle. 



and Vlll denotes projection of the wave function into 
the lowest Landau level. It was shown earlier that $1^ can 
be obtained from $„ by substituting the single electron 
wave functions Ya{vj) by the 'single CF wave functions' 

YS^{i^]) = 'Plll Y{k{zj-ZkYYo,{vj), the explicit form of 
which has been given in the literature pi|| . (The prime 
denotes the condition k ^ j.) The composite fermion 
interpretation of $1^ follows since multiplication by the 
Jastrow factor ni<fe(-^j ~ ZkY^ is tantamount to attach- 
ing 2p flux quanta to each electron, converting it into 
a composite fermion. These wave functions have been 
found to be quite accurate in tests against exact diago- 
nalization results available for small systems PJO]. 

The Hamiltonian for the many electron system is given 
by 



H=\Y.^{r,u) + Ve. 



(3) 
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where Ve-b is the electron-background interaction, with 
the background assumed to be comprised of a uniform 



positive charge, and V{r) is the effective two-dimensional 
electron-electron interaction. (The kinetic energy is 
quenched in the lowest Landau level.) For a strictly two- 
dimensional system, V{rj^.) — '^_ — r, where e is the 
dielectric constant of the background material. As we 
will see, an important quantitative correction comes from 
the finite transverse extent of the electron wave function, 
which alters the form of the effective two-dimensional 
interaction at short distances. The effective interaction 
can be calculated straightforwardly once the transverse 
wave function is known, which in turn will be deter- 
mined by self-consistently solving the Schrodinger and 
Poisson equations, taking into account the interaction 
effects through the local density approximation (LDA) 
including the exchange correlation potential |12| . Two 
geometries, single hetero junction and square quantum 
well (SQW), are considered due to their experimental rel- 
evance. To simplify the calculation, we assume that the 
electron wave function is confined entirely on the GaAs 
side of the heterojunction, which is a reasonably good 
approximation for deep confinement. It is stressed that 
neither the microscopic wave function nor the effective 
interaction contains any adjustable parameters; the for- 
mer depends only on the filling factor, while the latter 
is determined from a first principles, self-consistent LDA 
calculation, with the two-dimensional density, the sam- 
ple type (heterojunction or square quantum well) and the 
known sample parameters as the only input. 
The energy of the exciton sA v — 2n+i ' 



A'^ 



< $ff lif i$ff > < $rii?i*r > 






<$ri$r > 



(4) 



is computed by Monte Carlo methods in the standard 
spherical geometry [ [l3[ . Since moving a single particle 
at each step of the Monte Carlo changes the single CF 
wave function Y'^^ for all particles, the full wave func- 
tion must be computed at each step. The exciton wave 
function is a linear superposition of ~ N jn Slater deter- 
minants, but each of these differs from the ground state 
only in one row, and the clever techniques for upgrading 
Slater determinants [|l^ significantly reduce the comput- 
ing time, enabling us to study reasonably large systems 
(up to 63 composite fermions were used in the present 
study). The ground and excited state energies are evalu- 
ated sufficiently accurately to get a reasonable estimate 
for the gap. Due to the lack of edges in the spherical ge- 
ometry being studied, we expect that the gap will have 
a linear dependence on N^'^ to leading order, which is 
also borne out by our results. A linear extrapolation to 
the thermodynamic limit N^^ — » is taken after cor- 
recting the energies for the finite size deviation of the 
density from its thermodynamic value in the standard 
manner [|l5| . All results below are thermodynamic ex- 
trapolations, unless mentioned otherwise. The energies 
are quoted in units of e^ /d^, where Zq — \/he/ eB is the 



magnetic length. 

We have determined the dispersion of the CF exciton 
for 1/3, 2/5, and 3/7, corresponding to one, two, and 
three filled CF-LLs, respectively. The typical dispersion 
contains several minima, as shown in Fig. (nl). We will 
term the lowest energy minimum the "fundamental" ro- 
ton, or simply the roton, the others being secondary ro- 
tons. Since only discrete values of k are available in our 
finite systems, the energy of the roton is obtained by 
fitting the points near the minimum to a parabolic dis- 
persion 
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for each N, and then extrapolating A to the thermody- 
namic limit. The energies in the Mq — > limit and at 
the roton minimum are given in Table I for a strictly two 
dimensional system, along with m.^. 

Figs. (|^) and (^) plot the energy of the CF roton and 
the long wavelength CF exciton for a heterojunction and 
a square quantum well (of width 25nm) as a function 
of density, calculated with the realistic LDA interaction, 
along with the experimental energies obtained in phonon 
absorption (at 1/3, 2/5, and 3/7) as well as in inelastic 
light scattering experiments (at 1/3). A more detailed 
comparison is given in Table II. In the small wave vector 
limit, the calculated energy at 1/3 is off by ~ 30%. It 
has been suggested that here the true lowest energy exci- 
tation may contain two CF-excitons, and there has been 
debate as to which excitation is being probed by the Ra- 
man scattering in this case ||,|l^. For the roton, the 
theoretical energies, obtained with no adjustable param- 
eters, are in excellent agreement with the observed ones 
|17[ . One may worry that the situation will be spoiled 
by Landau level mixing. This turns out not to be the 
case. Following Ref. |jia], we have estimated the impor- 
tance of LL mixing for the roton energy by considering 
a variational wave function which is a linear combina- 
tion of the projected and unprojected wave functions, 
and found that the corrections are on the order of 5% for 
typical densities, consistent with similar conclusion for 
the transport gap [ p^ , p^ . 

While the ballistic phonon absorption experiments di- 
rectly measure the minimum energy (i.e., the roton) po| , 
the Raman experiments ideally probe the klo — > limit 
of the CF exciton dispersion, the wavelength of the light 
being much larger than Iq. However, a breakdown of mo- 
mentum conservation due to the presence of disorder can 
activate rotons as well |g^ , as a result of a singularity in 
the density of states. This has been crucial in explain- 
ing multiple peaks in the Raman spectra for the inter-LL 
excitations. At ^' = 1/3, a low energy Raman peak has 
been interpreted as the roton [g|4|. Recently, Kang et al. 
p2[ have also observed modes at 2/5 and 3/7, at energies 
of 0.031 and 0.008 e^/e^o, respectively, which they inter- 



pret as the long wavelength neutral mode. Our calcu- 
lated energies at 2/5 and 3/7 for a quantum well of width 
30nm and density p = 5.4 x lOi"^ cm'^ are 0.031(3) and 
0.021(3) e^/e^o, respectively, for the roton, and 0.070(1) 
and 0.056(3) e^/e/o for the klo = limit of CF exciton. 
At 2/5, the energy of the observed excitation is consis- 
tent not only with the calculated roton energy but also 
with those measured in ballistic phonon absorption ex- 
periments, and substantially smaller than the calculated 
Hq — > limit, which might suggest an identification with 
the roton. [We note here that the observation of the 1/3 
roton implies that the violation of the momentum con- 
servation is sufficiently wide-spread as to render the fun- 
damental rotons at 2/5 and 3/7 observable as well, which 
occur at roughly the same wave vectors {klo ~ 1-6 — 1.7) 
as the 1/3 roton {klo ~ 1-4).] The energy of the 3/7 
mode of Ref. p3] is anomalously low, however. Further 
work will be required to ascertain the origin of these new 
Raman modes; an experimental observation of multiple 
roton peaks will be especially helpful in clarifying this 
issue. 

In conclusion, the inscnsitivity of the roton energy to 
disorder has afforded an opportunity for a direct quan- 
titative confirmation of our theoretical understanding of 
the excitations of the fractional quantum Hall effect. We 
are grateful to Professor Aron Pinczuk for communicat- 
ing his results to us prior to publication and for the con- 
tinuous exchange of information, and to Xiaomin Zu for 
numerous helpful discussions. This work was supported 
in part by the National Science Foundation under grant 
no. DMR-9615005, and by a grant of computing time by 
the National Center for Supcrcomputing Applications at 
the University of Illinois (Origin 2000). 
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TABLE I. Energies of the CF roton and the long wave- 
length neutral excitation at 1/3, 2/5, and 3/7, for a strictly 
two-dimensional system, in units of e^ /do. Also given is an es- 
timate for the CF roton "mass", m|j = /i|jm,e\/-B[T], defined 
in Eq. (pi), where me is the electron mass in vacuum. The sta- 
tistical uncertainty in the last digit(s) is shown in parentheses. 



mode 


V 


energy 


Mh 


klo =0 


1/3 


0.15 


- 


2/5 


0.087(1) 


- 


3/7 


0.068(5) 


- 


roton 


1/3 


0.066(1) 


0.0079(3) 


2/5 


0.037(1) 


0.0090(11) 


^ 


0.027(3) 


0.0095(32) 




TABLE II. Comparison of theory and experiment for the 
roton energy as well as the energy of the long- wavelength ex- 
citon, quoted in units of e^ jdo- In Ref. H, the roton energies 
were determined for 2/3, 3/5, and 4/7, which, assuming parti- 
cle-hole symmetry, are the same as the roton energies at 1/3, 
2/5, and 3/7, when measured in units of e^/e/o- 
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klo = 





roton 


Reference 


experiment 


theory 


experiment 


theory 






1/3 


0.082 


0.104(1) 


0.044 


0.050 (1) 






0.084 


0.113(1) 


- 


0.052(1) 






- 


0.09(2) 


0.041(2) 


0.045(1) 






0.074 


0.095(1) 


0.047 


0.047(1) 






- 


0.092(1) 


0.036(5) 


0.045(1) 






2/5 


- 


0.054(1) 


0.021(2) 


0.026(1) 






- 


0.055(1) 


0.025(3) 


0.027(1) 






3/7 


- 


0.044(2) 


0.014(2) 


0.017(2) 







FIG. 1. The dispersions of the CF excitons at // = 1/3, 2/5, 
and V = 3/7 for a strictly two dimensional system. The curves 
are obtained from the discrete points of systems with up to 50 
particles (without extrapolation to the thermodynamic limit), 
with the typical Monte Carlo uncertainty shown at the begin- 
ning of each curve. The energies are given in units of e^/eZo, 
where e = 12.8 is the dielectric constant of GaAs, and lo is 
the magnetic length. 
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FIG. 2. The energies of the CF roton (sohd line) and the 
long-wavelength CF exciton (dashed line) for heterojunction 
as a function of density, p. The typical uncertainty in the 
theoretical results is shown at the beginning of each curve. 
Experimental energies are also shown, taken from Ref. [i] 
(circle), Ref. M (diamond), Ref. H (square), and Ref. [?] 
(down-triangle); the filled symbols correspond to the roton, 
and the empty ones to the long wavelength mode. 



FIG. 3. Same as in Fig. (Q) but for square quantum well of 
width 25nm. The experimental results are taken from Ref. M 
(up-triangle) and Ref. M (right-triangle); the filled (empty) 
symbols correspond to the roton (long wavelength mode). 



